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Preface 


This  report  illustrates  the  evolution  of  directionally  distributed,  shallow- 
water  waves  at  times  of  high  energy  as  measured  with  a  high-resolution  direc¬ 
tional  wave  gauge.  The  work  was  motivated  by  the  need  for  a  complete  de¬ 
scription  of  storm  wave  fields  incident  on  the  nearshore  to  augment  under¬ 
standing  and  modeling  of  nearshore  processes  which  affect  coastal  engineering 
projects.  This  effort  was  authorized  by  Headquarters,  U.S.  Army  Corps  of 
Engineers  (HQUSACE),  under  Civil  Works  Research  Work  Unit  32484, 
Directionality  of  Waves  in  Shallow  Water,  Coastal  Flooding  Program.  Funds 
were  provided  through  the  Coastal  Engineering  Research  Center  (CERC), 

U.S.  Army  Engineer  Waterways  Experiment  Station  (WES),  under  the  pro¬ 
gram  management  of  Ms.  Carolyn  M.  Holmes,  CERC.  Messrs.  John  H. 
Lockhart,  Jr.,  John  G.  Housley,  Barry  W.  Holliday,  and  David  A.  Roellig 
were  HQUSACE  Technical  Monitors. 

This  report  was  prepared  by  Dr.  Charles  E.  Long  at  CERC’s  Field  Re¬ 
search  Facility  (FRF)  in  Duck,  NC,  under  the  direct  supervision  of 
Mr.  William  A.  Birkemeier,  Chief,  FRF,  and  Mr.  Thomas  W.  Richardson, 
Chief,  Engineering  Development  Division,  CERC;  and  under  the  general 
supervision  of  Dr.  James  R.  Houston  and  Mr.  Charles  C.  Calhoun,  Jr.,  Di¬ 
rector  and  Assistant  Director,  CERC,  respectively. 

Data  processing  software  was  written  by  Dr.  Joan  M.  Oltman-Shay  of 
Quest  Integrated,  Inc.,  Seattle,  WA.  Physical  maintenance  of  the  directional 
gauge  was  done  by  the  FRF  dive  team  consisting  of  Messrs.  Birkemeier,  Mi¬ 
chael  W.  Leffler,  H.  Carl  Miller,  Eugene  W.  Bichner,  and  Brian  L. 
Scarborough.  Gauge  calibration  was  maintained  by  Mr.  Kent  K.  Hathaway, 
FRF.  Data  acquisition  and  storage  were  done  by  Mr.  Clifford  F.  Baron, 

FRF.  Considerable  assistance  in  data  processing  was  provided  by  students 
contracted  through  the  Cooperative  Education  Program:  Ms.  Wendy  L.  Smith 
of  Old  Dominion  University,  Mmes.  Juliana  Admadji  and  Janna  Pemberton  of 
Florida  Institute  of  Technology,  and  Ms.  Judy  Roughton  of  College  of  the 
Albemarle,  Elizabeth  City,  NC.  Dr.  R.  T.  Guza,  Scripps  Institution  of 
Oceanography,  and  Mr.  Lockhart  provided  several  useful  suggestions  after 
reviewing  a  preliminary  draft  of  this  manuscript.  The  efforts  of  all  of  these 
individuals  are  appreciated. 

At  the  time  of  publication  of  this  report,  Director  of  WES  was  Dr.  Robert 
W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


1  Introduction 


Ocean  wind  waves  are  among  the  dominant  forcing  mechanisms  for  near¬ 
shore  dynamic  processes.  Of  the  broad  suite  of  possible  wave  regimes,  those 
associated  with  high  energy  are  most  important  for  engineering  design.  High- 
energy  waves  produce  the  most  rapid  changes  in  natural  or  renourished  beach 
morphologies  and  are  most  likely,  by  definition,  to  exceed  the  strength  or 
elevation  of  a  coastal  protective  structure,  resulting  in  costly  damage  or  exces¬ 
sive  overtopping.  Complete  knowledge  of  the  types  of  extremely  energetic 
natural  wind  wave  fields  that  are  incident  on  the  nearshore  is  thus  rather 
critical  for  efficient  design  of  coastal  engineering  projects. 

A  complete  wave  field  definition  describes  the  way  in  which  wave  energy 
is  distributed  in  both  frequency  and  direction.  Because  of  the  difficulty  in 
making  high-resolution  measurements  of  wave  energy  distributions,  early 
engineering  guidance  treated  a  given  sea  state  as  consisting  of  a  train  of  waves 
at  a  single  frequency  traveling  in  a  single  direction.  While  most  wave  fields 
are  more  complicated  than  this,  there  persists  much  guidance  based  on  unidi¬ 
rectional,  monochromatic  waves,  as  evidenced  in  the  Shore  Protection  Manual 
(SPM  1984).  Advances  in  time  series  measurement  and  efficient  analysis  in 
the  last  three  decades  have  provided  considerable  knowledge  of  the  frequency 
distributions  of  wave  energy,  enabling  the  development  of  good  theoretical 
models  and  improved  physical  models  of  wave  processes  with  a  consequent 
improvement  in  engineering  guidance. 

Detailed  measurements  of  wave  direction  are  more  difficult.  Only  within 
the  last  decade  have  practical  methods  of  directional  wave  measurement  been 
available,  and  these  are  still  expensive  and  logistically  demanding.  Until  very 
recently,  there  were  no  long-term  observations  from  which  to  extract  the 
directional  character  of  nearshore  seas  with  any  statistical  confidence.  To 
alleviate  this  lack  of  knowledge,  the  Field  Research  Facility  (FRF)  of  the  U.S. 
Army  Engineer  Waterways  Experiment  Station,  Coastal  Engineering  Research 
Center,  deployed  a  high-resolution  directional  wave  gauge  at  its  site  near 
Duck,  NC,  in  September  1986.  Measurements  covering  a  5-year  period, 
lasting  through  August  1991,  enable  considerable  clarification  of  the  direc¬ 
tional  nature  of  incident  wind  waves  at  this  site. 

The  purpose  of  this  report  is  to  elucidate,  by  illustration,  the  nature  of  the 
frequency  and  directional  distributions  of  wind  wave  energy  during  storm 
events  occurring  during  the  5-year  observation  period.  These  results  are 


Chapter  1  Introduction 


intended  to  assist  numerical  modelers  in  developing  and  testing  theories  of 
wave  transformation  across  a  continental  shelf.  Results  also  can  be  used  by 
numerical  and  physical  modelers  to  establish  realistic  seaward  wind  wave 
boundary  conditions  for  models  of  nearshore  processes.  Illustrations  are  con¬ 
tained  primarily  in  Appendix  A.  Chapters  2  and  3  of  this  report  describe  the 
experiment  site,  gauges,  sampling  pattern,  storm  definition,  characterizing 
parameters,  and  for  nat  of  displayed  results.  Common  features  of  storm  wave 
frequency-direction  spectra  are  described,  and  some  conjecture  about  unusual 
spectral  properties,  along  with  the  appropriate  path  of  future  research,  is 
presented  in  Chapter  4. 
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2  Measurement  Scheme 


Measurement  Site 

The  FRF  is  located  on  the  Outer  Banks  of  North  Carolina,  just  north  of  the 
village  of  Duck  (Figure  1).  The  site  is  representative  of  many  coastal  sites 
characterized  by  a  broad  continental  shelf,  relatively  mild  bottom  slope,  and 
reasonably  uniform  longshore  topography.  The  mild  sloping  bottom  substan¬ 
tially  dissipates  wind  waves  in  the  nearshore  so  that  attention  can  be  focused 
on  waves  arriving  from  seaward.  The  uniform  longshore  topography  allows 
the  use  of  a  spatial  array  of  sensors  to  detect  wave  direction  without  contami¬ 
nation  by  a  longshore  inhomogeneity  in  the  incident  wave  field.  The  facility, 
described  by  Birkemeier  et  al.  (1985),  has  a  permanent  staff  and  computing 
facilities  that  are  used  to  monitor,  maintain,  and  error  check  the  instrumenta¬ 
tion  on  a  day-to-day  basis,  and  so  minimize  data  losses  from  the  myriad  prob¬ 
lems  that  can  plague  field  data  collection. 

Regional  and  local  topography 

At  the  large  environmental  scale,  the  coastline  in  the  vicinity  of  the  FRF  is 
nearly  straight  for  several  tens  of  kilometers  north  and  south.  It  is  oriented 
such  that  the  shore-normal  direction  (facing  seaward)  is  very  nearly  70  deg 
from  true  north  (70  deg  T).  Waves  and  onshore  winds  can  approach  this  site 
along  an  easterly  180-deg  arc  from  340  to  160  deg  T.  The  adjacent  continen¬ 
tal  shelf  is  approximately  100  km  wide,  being  narrower  to  the  south  where 
Cape  Hatteras  juts  eastward.  The  direction  of  nearest  approach  of  the  100-m 
isobath  (the  shelf  break)  is  10  deg  to  15  deg  south  of  east  and  is  about  80  km 
distant  from  the  FRF.  A  characteristic  bottom  slope  for  the  shelf  is  1  m/km, 
but,  at  a  finer  scale,  the  bathymetry  is  marked  by  numerous  features  of  1-  to 
10-km  horizontal  scales  and  10-m  vertical  scales  scattered  irregularly  across 
the  shelf.  Waves  with  periods  near  10  sec  begin  to  shoal  at  depths  of  100  m, 
so  that  considerable  refraction  of  wind  waves  propagating  over  the  shelf  is 
expected  (U.S.  Army  Engineer  District,  Wilmington  1980). 

Within  a  few  kilometers  of  the  FRF,  the  bathymetry  is  more  regular.  A 
complex  bar  system  exists  within  about  300  m  of  the  shore  (Birkemeier  1984), 
and  waves  and  currents  have  created  some  irregular  bathymetry  of  the  600-m- 
long  research  pier  (Miller,  Birkemeier,  and  DeWall  1983).  Away  from  these 
regions,  isobaths  are  nearly  shore-parallel.  The  bottom  slope  is  nearly  uni¬ 
form  at  about  0.5  m/km  and  is  reasonably  dynamically  stable,  as  indicated  by 
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Figure  1 .  FRF  location  and  offshore  bathymetry 

repeated  surveys  (Howd  and  Birkemeier  1987).  Figure  2  illustrates  the 
bathymetry  over  a  region  600  m  north  to  600  m  south  of  the  pier  and  extend¬ 
ing  about  900  m  offshore. 

Regional  climate 

The  site  is  subject  to  a  variety  of  climates  that  give  rise  to  a  broad  diversi¬ 
ty  in  frequency  and  directionally  distributed  wave  energy.  Typical  wind  and 
wave  climate  for  the  period  covered  by  this  report  can  be  roughly  classified  as 
four  basic  types  (Miller  et  al.  1988;  Leffler  et  al.  1989,  1990,  1991,  1992, 
1993).  First  are  low-wind  situations,  where  the  wave  field  is  characteristically 
low-frequency  swell  propagating  in  from  the  open  Atlantic  Ocean.  These 
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Figure  2.  FRF  nearshore  bathymetry  and  instrument  location 


waves  typically  arrive  from  the  east-to-southeast,  and  are  normally  present 
between  storm  events. 

The  second  and  most  common  type  of  wave-generating  winds  are  from 
northeasters  associated  with  atmospheric  frontal  passages.  Fronts  typically 
arrive  from  the  west  and  have  a  northeast-southwest  orientation.  Active  local 
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wave  generation  occurs  after  a  frontal  passage  when  strong  winds  arise  in  the 
northeast  quadrant.  A  third  type  of  storm  is  an  extratropica)  depression  that 
typically  forms  south  of  the  FRF  along  the  south  Atlantic  coast.  These  storms 
tend  to  migrate  northward  along  the  coast,  creating  high  winds  and  waves 
from  the  east  to  northeast.  Les:  frequent  but  more  intense  wave  generators 
are  hurricanes,  the  fourth  type  of  disturbance.  These  originate  in  the  equatori¬ 
al  Atlantic,  migrate  west  to  northwest,  and  can  make  landfall  or  the  gulf  and 
Atlantic  coasts  with  very  destructive  winds  and  waves.  No  hurricanes  made 
landfall  at  the  FRF  during  the  course  of  this  study,  but  several  passed  by  in 
the  deeper  Atlantic,  creating  a  distinctive  local  field  of  high,  long-period 
waves  narrowly  distributed  in  frequency  and  direction. 

It  is  noted  that  the  regional  topography  and  climate  limit  the  generalization 
of  results  obtained  at  the  FRF  to  sites  having  similar  characteristics.  For 
example,  one  would  not  expect  nearshore  dynamics  of  sites  on  the  Great 
Lakes,  gulf,  or  west  coasts  of  the  United  States  to  behave  in  the  same  way  be¬ 
cause  their  climates  and  bathymetries  are  different.  Nonetheless,  fundamental 
physical  processes  such  as  nonlinear  wave  behavior,  wave  breaking,  and  air- 
sea  interaction,  for  which  insights  are  gained  by  studies  at  the  FRF,  will  be 
more  universally  applicable;  the  limitations  are  not  too  restrictive,  and  gener¬ 
alizations  are  justified. 

Instrumentation 

Data  from  two  gauges  are  presented  in  this  report.  One  gauge  consists  of 
a  set  of  near-bottom-mounted  pressure  gauges,  which  together  constitute  a 
high-resolution  directional  wave  gauge.  This  gauge  provides  estimates  of  the 
frequency-direction  spectra  of  wind  waves  which  are  the  basis  of  this  report. 
The  other  gauge  is  an  anemometer  that  gives  a  measure  of  local  wind  speed 
and  direction.  It  provides  a  characterization  of  the  larger  wind  field  that 
generates  some  of  the  wave  fields  measured  with  the  directional  wave  gauge. 
These  gauges  are  described  here  briefly  for  completeness  and  because  their 
properties  affect  the  way  data  are  interpreted. 

High-resolution  directional  wave  gauge 

The  directional  wave  gauge  consists  of  two  fundamental  parts:  (a)  an  array 
of  gauges  that  sense  sea  surface  displacement,  and  (b)  an  analysis  algorithm 
that  converts  the  data  to  an  estimate  of  S(f,6),  the  frequency-direction  spec¬ 
trum,  where  /  is  cyclic  wave  frequency  and  6  is  wave  direction.1  The  sen¬ 
sors  are  a  set  of  nine  near-bottom-mounted  pressure  gauges,  laid  out  in  a 
straight  line  along  the  8-m  isobath  as  shown  in  Figure  2.  Spacing  between  the 
gauges  follows  the  guidance  given  by  Davis  and  Regier  (1977)  for  the  design 
of  wave  direction  detectors.  Pressure  time  series  data  from  each  gauge  are 
converted  to  units  of  equivalent  seawater  depth,  Fourier  transformed,  and 
converted  to  sea  surface  displacement  Fourier  transforms  using  linear  wave 


1  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation  (Appendix  B). 
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theory.  For  each  pair  of  gauges,  raw  cross-spectral  estimates  are  then  com¬ 
puted,  ensemble  averaged,  and  band  averaged  to  create  a  smoothed  cross¬ 
spectrum  for  each  gauge  pair.  For  each  smoothed  frequency,  a  matrix  is 
formed  of  the  cross-spectra  for  all  possible  pairs  of  gauges.  This  matrix  is 
then  fed  to  the  directional  estimation  algorithm. 

Wave  energy  directional  distributions  are  computed  using  the  Iterative 
Maximum  Likelihood  Estimator  (IMLE),  which  was  derived  by  Pawka  (1982, 
1983),  and  found  to  be  a  substantial  improvement  on  the  Maximum  Likeli¬ 
hood  Estimator  (MLE)  derived  by  Davis  and  Regier  (1977).  Specifically,  the 
algorithm  is  Equations  1  -  4  in  the  paper  by  Pawka  (1983)  with  convergence 
exponent  jS  =  1 ,  convergence  coefficient  y  =  5 ,  and  a  maximum  of  30  itera¬ 
tions. 

In  practice,  all  gauges  of  the  FRF  array  are  not  used  for  all  wave  frequen¬ 
cies.  Gauges  with  large  separations  are  used  for  long,  low-frequency  waves. 
Closely  spaced  gauges  are  used  for  shorter  waves  to  avoid  loss  of  coherence 
for  these  waves  at  large  spatial  separations.  Details  of  array  sub-configura¬ 
tions  and  their  associated  frequencies  are  given  by  Long  and  Oltman-Shay 
(1991).  Tests  of  the  FRF  array  with  artificial  data  indicate  a  directional  reso¬ 
lution  of  5  deg.  This  means  it  can  resolve  the  peak  of  a  single  mode  direc¬ 
tional  distribution  to  within  S  deg,  or  distinguish  the  peaks  of  a  bimodal  distri¬ 
bution  when  separated  by  15  deg  or  more. 

Anemometer 

The  reference  anemometer  used  in  this  study  was  mounted  atop  the  FRF 
laboratory  building  (Figure  2)  at  an  elevation  of  19.5  m  above  mean  sea  level. 
It  has  a  nominal  speed  accuracy  of  ±0.5  m/sec  and  direction  accuracy  of 
±3  deg.  Speed  and  direction  data  are  collected  on  separate  channels,  and 
results  from  each  channel  are  simply  averaged  (not  vector  averaged).  The 
proximity  of  the  anemometer  to  the  FRF  building  is  known  to  cause  low  speed 
readings  on  some  occasions  in  very  high  winds  (Long  and  Hubertz  1988),  but 
it  is  considered  sufficiently  reliable  for  the  generally  moderate  winds  encoun¬ 
tered  in  this  study.  Though  other  anemometers  have  been  deployed  for  short¬ 
er  periods  at  the  FRF,  this  one  was  consistently  in  use  for  the  duration  of  this 
study.  It  follows  that  this  gauge  provides  the  most  consistent  local  represen¬ 
tation  of  the  larger  scale  wind  field. 

Data  Sampling  Schemes 

With  few  exceptions,  data  used  in  this  report  were  obtained  following  the 
standard  FRF  sampling  scheme.  The  pattern  is  to  collect  four  contiguous 
records  of  34-min  duration  (4,096  data  points  at  a  2-Hz  sampling  rate)  starting 
at  0100,  0700,  1300,  and  1900  Eastern  Standard  Time  (EST).  When  the 
characteristic  wave  height  exceeds  2  m  at  a  reference  wave  gauge  (usual¬ 
ly  gauge  tt\  of  the  directional  array),  additional  collections  are  made  starting 
at  0400,  1000,  1600,  and  2200  EST.  This  form  of  sampling  gives  good  data 
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coverage  of  high  energy  conditions,  but  has  a  tendency  to  miss  some  of  the 
details  of  initial  transient  conditions  at  the  onset  of  a  storm. 

The  anemometer  data  presented  in  this  report  are  based  on  a  subset  of  the 
overall  FRF  database.  Wind  speeds  and  directions  are  deduced  from  single 
34-min  records  with  start  times  of  0100,  0700,  1300,  and  1900  EST.  These 
particular  records  are  the  most  reliable,  having  been  checked  by  hand  for  use 
in  the  FRF  annual  reports  (Leffler  et  al.  1992,  1993).  The  result  is  a  very 
regular  sampling  of  local  winds,  with  gaps  only  for  system  maintenance  or 
repair. 

All  possible  data  were  processed  from  the  directional  wave  gauge,  so  that 
normally  there  are  samples  from  all  the  routine  start  times  listed  above.  Four 
records  of  34  min  each  (2  hr  16  min  or  16,384  data  points  at  2  Hz)  for  each 
gauge  were  used  in  cross-spectral  computation  to  ensure  high  confidence  in 
the  resulting  cross-spectral  matrix  that  is  transformed  in  the  IMLE  algorithm. 
Fifteen  half-overlapping  segments  of  2,048  points  were  processed  for  ensem¬ 
ble  averages,  and  10  subsequent  frequency  bands  were  averaged  for  a  final 
resolution  bandwidth  A/  =  0.00977  Hz  and  a  minimum  of  160  degrees  of 
freedom  (excluding  any  gain  in  degrees  of  freedom  from  the  lapped  segments) 
in  the  final  cross-spectral  estimates.  While  this  method  gives  very  good  re¬ 
sults  in  estimates  of  S(f,6) ,  it  has  the  effect  of  averaging  all  wave  processes 
for  the  2  hr  16  min  records  so  that  highly  transient  processes  are  smoothed 
somewhat.  Results  are  not  invalid,  just  rather  strongly  filtered.  Shorter 
records  can  be  used  in  future  work  if  brief  transients  are  important. 

The  bulk  of  S{f,6 )  results  were  obtained  by  the  method  just  described,  but 
there  were,  on  occasion,  minor  variations  in  sampling  times  and  rates.  Dur¬ 
ing  the  SUPERDUCK  experiment  in  September  and  October  1986  (Crowson 
et  al.  1988),  data  collection  start  times  were  adjusted  to  conform  to  high  and 
low  tides,  and  so  did  not  conform  to  the  FRF  normal  collection  pattern. 
Furthermore,  a  total  of  seven  34-min  records  were  obtained  in  each  collection. 
To  maximize  the  number  of  S(f,6)  observations  during  SUPERDUCK,  re¬ 
cords  1  -  4  from  a  collection  were  processed,  and  then  records  4  -  7  were 
processed,  so  that  record  4  was  used  twice,  but  observations  were  still  almost 
independent.  Time  series  of  parameters  from  these  spectra  thus  appear  rather 
oddly  spaced,  with  pairs  of  samples  appearing  with  1  hr  42  min  separations. 
During  the  DELILAH  experiment  of  October  1990  (Birkemeier  et  al.  1991), 
data  were  collected  for  an  additional  34  min  and  the  sampling  rate  for  all  data 
was  4  Hz  instead  of  2  Hz.  Thus,  a  collection  had  10  records  of  4,096  data 
points  (about  17  min)  each.  For  the  data  used  here,  wind  speeds  and  direc¬ 
tions  during  DELILAH  were  averaged  for  only  17  min  instead  of  the  usual 
34  min.  Directional  array  data  were  still  processed  in  2  hr  16  min  total  re¬ 
cord  lengths,  but,  in  smoothing,  20  bands  were  averaged  instead  of  10,  result¬ 
ing  in  the  same  final  resolution  bandwidth  (but  more  degrees  of  freedom). 
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General  Working  Data  Set 


The  final  working  data  set  from  which  to  extract  storm  data  consists  of 
6,762  observations  of  S(f,6)  spanning  the  period  from  September  1986  to 
August  1991,  or  about  5  years.  Two  major  gaps  in  data  collection  occurred 
during  this  time,  affecting  the  seasonal  coverage  of  storms  and  reducing  their 
total  number.  After  SUPERDUCK,  full  time  collection  of  directional  gauge 
data  did  not  resume  until  late  January  1987,  so  there  is  a  gap  in  the  time  line 
of  collections  from  November  1986  to  January  1987.  In  December  1987,  a 
fishing  boat  trawled  a  net  along  the  line  of  the  linear  array  of  gauges,  destroy¬ 
ing  about  half  of  them.  Parts  replacement  and  repairs  could  not  be  effected 
until  August  1988,  so  there  is  another  gap  from  December  1987  to  August 
1988.  Aside  from  these  two  major  gaps,  the  record  is  substantially  complete, 
excepting  short  periods  for  routine  maintenance  and  calibration. 
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3  Storm  Data 


Event  Definition 


To  extract  storm  data  from  the  larger  working  data  set,  a  definition  of  what 
is  meant  by  a  storm  is  necessary.  For  the  purposes  of  this  report,  a  storm  is 
defined  using  the  characteristic  wave  height  and  the  length  of  time  it 
exceeds  a  certain  level.  To  isolate  the  largest  events,  it  was  required  that 
exceed  2  m  for  times  at  least  24  hr  apart,  and  not  drop  below  1  m  in  the 
interim.  Wave  conditions  are  normally  relatively  quiescent  just  before  a 
storm,  generally  being  low-amplitude,  long-period  swell  from  the  southeast 
quadrant.  To  ensure  at  least  one  prestorm  condition,  data  for  a  given  event 
were  extracted  from  the  main  database  from  the  beginning  of  the  day  contain¬ 
ing  a  prestorm  condition  to  the  end  of  the  day  during  which  first  drops 
below  1  m. 

Figure  3  shows  a  time  series  of  H ^  from  the  directional  gauge  for  all  col¬ 
lections  in  the  main  database.  Storm  events  that  satisfy  the  above  definition 
are  identified  by  letter  or  number.  Their  nominal  durations  are  marked  by 
horizontal  black  bars.  In  all,  29  storm  events  were  found  and  are  referenced 
by  the  letters  A  to  Z  and  the  numbers  1  to  3.  Note  that  Event  3  follows 
Event  E. 

Table  1  lists  the  storm  events  shown  in  Figure  3,  along  with  the  start  and 
end  dates  defined  by  the  criteria  described  above.  Also  shown  are  the  maxi¬ 
mum  H ^  observed  during  the  storm  and  a  nominal  duration,  defined  as  the 
time  elapsed  from  the  first  time  H ^  exceeds  2  m  to  the  last  time  it  drops 
below  2  m.  Linear  interpolation  between  measured  H ^  was  used  to  obtain 
the  exceedence  times.  Based  on  the  definition  of  a  storm  described  above,  all 
durations  should  exceed  24  hr.  Several  do  not,  which  illustrates  the  pitfalls  of 
searching  a  6,762-point  database  by  eye.  The  shorter  storms  are  retained  here 
because  they  do  not  differ  from  the  longer  storms  except  by  their  durations. 

Climatological  Parameters 

Several  characterizing  parameters  are  used  in  this  report  to  describe  a 
storm  evolution.  These  parameters  are  given  names  and  symbols  common  to 
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Figure  3.  Time  series  of  Hmo  with  storm  events  identified 
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Table  1 

Storm  Event  Identification 

Event 

Start 

End 

m 

duration 

hr 

A 

9 

Oct 

86 

14 

Oct 

86 

3.09 

57.01 

B 

17 

Oct 

86 

21 

Oct 

86 

2.38 

23.67 

C 

15 

Feb 

87 

20 

Feb 

87 

4.03 

40.07 

D 

9 

Mar 

87 

18 

Mar 

87 

4.22 

161.84 

E 

22 

Mar 

87 

26 

Mar 

87 

2.68 

28.15 

3 

11 

Oct 

87 

17 

Oct 

87 

3.11 

85.29 

F 

10 

Nov 

87 

12 

Nov 

87 

2.45 

27.04 

G 

26 

Nov 

87 

1 

Dec 

87 

2.21 

86.13 

H 

14 

Dec 

88 

17 

Dec 

88 

2.31 

27.34 

1 

3 

Jan 

89 

6 

Jan 

89 

2.19 

8.68 

J 

20 

Jan 

89 

27 

Jan 

89 

3.00 

78.70 

K 

16 

Feb 

89 

20 

Feb 

89 

3.19 

54.43 

L 

22 

Feb 

89 

26 

Feb 

89 

4.32 

52.32 

M 

2 

Mar 

89 

5 

Mar 

89 

2.11 

24.45 

N 

6 

Mar 

89 

14 

Mar 

89 

4.23 

108.62 

O 

21 

Mar 

89 

25 

Mar 

89 

2.53 

30.08 

P 

3 

Sep 

89 

11 

Sep 

89 

2.59 

138.10 

Q 

18 

Sep 

89 

28 

Sep 

89 

2.68 

199.53 

R 

22 

Oct 

89 

1 

Nov 

89 

2.63 

38.02 

S 

7 

Dec 

89 

11 

Dec 

89 

4.17 

61.76 

T 

21 

Dec 

89 

26 

Dec 

89 

4.68 

91.00 

U 

6 

Mar 

90 

9 

Mar 

90 

2.96 

22.25 

V 

25 

Oct 

90 

28 

Oct 

90 

4.70 

44.31 

w 

17 

Nov 

90 

21 

Nov 

90 

2.51 

42.70 

X 

7 

Dec 

90 

10 

Dec 

90 

2.41 

24.07 

Y 

7 

Jan 

91 

14 

Jan 

91 

3.40 

111.84 

z 

18 

Apr 

91 

21 

Apr 

91 

3.38 

51.00 

1 

18 

May 

91 

21 

May 

91 

2.82 

33.22 

2 

23 

Jun 

91 

28  Jun 

91 

2.66 

21.05 

wave  analysis,  but  because  there  are  several  ways  to  compute  each  with  the 
data  at  hand,  they  are  defined  specifically  here.  The  parameter  set  consists  of 
characteristic  wave  height  H^,  characteristic  wave  period  Tp,  characteristic 
wave  direction  6p,  representative  directional  spread  of  wave  energy  86 p,  and 
water  depth  d. 

Characteristic  wave  height  H ^  is  defined,  as  usual,  as  four  times  the 
standard  deviation  of  sea  surface  displacement,  based  on  variances  contained 
in  the  frequency  spectrum  S(f)  between  the  low-  and  high-frequency  cutoffs. 
In  this  analysis,  S(J)  is  computed  as  the  average  of  the  surface-corrected 
frequency  spectra  from  each  of  the  nine  gauges  in  the  linear  array.  This  does 
not  increase  degrees  of  freedom  (because  all  gauges  are  measuring  the  same 
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wave  field),  but  it  does  reduce  instrument  noise.  All  analysis  is  done  in  dis¬ 
crete  frequency  space,  and  results  are  defined  in  terms  of  N  *  28  frequency 
bands  from/,  =  0.054  Hz  (the  low-frequency  bound)  to  fN  =  0.319  Hz  (the 
high-frequency  bound)  in  steps  of  A/  =  0.00977  Hz.  The  discrete  frequency 
spectrum  is  given  the  symbol  S(fn)  so  that  characteristic  wave  height  is 


-  4 


E 


to 


Characteristic  wave  period,  or  peak  period  Tp  is  defined  as  the  inverse  of 
the  peak  frequency^  deduced  from  the  frequency  spectrum;  i.e.,  the  discrete 
frequency  at  which  S{fJ  is  maximum.  This  is  the  conventional  definition  of 
peak  period,  and  is  used  here  to  be  consistent  with  other  wave  studies.  Thus, 
Tp  -  1  lfp. 

There  are  a  number  of  ways  to  define  peak  direction  6p.  The  one  used 
here  is  the  peak  of  the  discrete  frequency-direction  spectrum  S(fn,6J  at  the 
peak  frequency;  i.e.,  the  discrete  direction  6p  at  which  S(fp,6m)  is  a  maxi¬ 
mum.  Discrete  directions  6m  range  from  0,  =  -90  deg  to  6U  =  09l  =  90  deg 
in  A0  =  2-deg  increments.  These  angles  represent  the  directions  from  which 
wave  energy  arrives  relative  to  shore-normal,  with  -90  deg  being  waves  trav¬ 
elling  parallel  to  the  beach  from  the  south,  0  deg  being  waves  arriving  along 
the  shore-normal  axis,  and  +90  deg  being  waves  travelling  parallel  to  the 
beach  from  the  north.  A  constraint  on  linear  array  analysis  is  that  it  cannot 
detect  waves  travelling  seaward,  so  all  waves  are  assumed  to  arrive  from  the 
seaward  semicircle  of  azimuths. 


A  characteristic  directional  spread  of  wave  energy  5 6p  also  is  a  very  useful 
parameter  for  these  data.  It  is  also  defined  from  S(fn,6J  as  being  the  angle 
that  subtends  the  central  half  of  the  energy  in  the  directional  distribution  at  the 
peak  frequency.  Symbolically,  an  integral  function  l(f,6)  can  be  defined  as 


m  =  | 


6 

-90 


sm 

SV) 


d<t> 


(2) 


where  Iff,- 90°)  =  0  and  Itf,90°)  =  1  identically  because  the  integral  of 
S(f>6)  with  respect  to  direction  defines  S(f) .  The  arc  that  contains  the  central 
half  of  the  energy  at  a  given  frequency  begins  at  the  angle  0^  for  which 
l(f,6 jj)  =  0.25  and  ends  at  the  angle  075  for  which  l(f,6ls)  =  0.75.  If  0^ 
and  07J  are  determined  at  /  =  fp,  the  peak  directional  spread  is  simply 


60,  =  07J  -  0* 

This  parameter  is  robust  in  that  it  can  be  computed  for  any  finite  directional 
distribution  and  is  meaningful  in  that  it  gives  a  small  number  for  narrowly 
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distributed  energy,  and  a  large  number  for  broadly  distributed  energy.  Note 
that  its  use  here  differs  from  that  of  Long  and  Oltman-Shay  (1991).  Long  and 
Oltman-Shay  (1991)  looked  at  the  overall  spread  of  Examined  here  is 

the  energy  spread  only  at  /  ■  f ,  the  most  energetic  frequency  spectral  line. 

Water  depths  d  reported  here,  and  used  for  tide  stage  determination  and 
surface  correction  of  near-bottom  pressure  signals,  are  the  medians  of  the 
mean  depths  determined  from  the  nine  pressure  gauges  for  the  duration  of  one 
collection.  This  was  necessary  because,  while  the  pressure  gauges  are  accu¬ 
rate  at  wind  wave  frequencies,  they  are  susceptible  to  long-term  drifts  in  their 
mean  signals  due  to  environmental  changes  and  biological  activity  within  the 
gauge  sensing  volumes.  Use  of  median  depths  as  representative  depths  for  all 
the  gauges  eliminates  the  largest  deviate  mean  depths.  Comparison  of  result¬ 
ing  mean  depths  with  the  FRF  tide  gauge  indicates  an  accuracy  of  about 
0. 1  m,  which  is  sufficient  for  use  here. 

Data  Presentation 

The  primary  result  of  this  study  is  the  set  of  figures  contained  in  Appendix 
A.  For  each  of  the  events  listed  in  Table  1,  there  are  two  types  of  display. 
One  page  contains  time  series  of  climatological  parameters  associated  with  the 
event.  Following  the  climatology  page  for  each  event,  subsequent  pages 
contain  sequences  of  frequency-direction  spectral  plots  for  select  times  within 
the  event.  Page  A3  contains  an  index  of  events,  showing  the  page  in  Appen¬ 
dix  A  where  climatological  parameters  are  displayed. 

Climatological  parameters  displayed  for  the  duration  of  each  event  are 
those  described  above,  namely  H^,  Tp,  d,  5 6p,  6p,  wind  direction,  and  wind 
speed.  The  reader  is  referred  to  the  parametric  plot  for  Event  A  on  page  A4 
for  an  example.  Wind  direction  and  peak  wave  direction  8p  have  been  trans¬ 
formed  to  the  same  coordinate  system  (directions  from  which  winds  or  waves 
arrive  in  degrees  counterclockwise  from  shore-normal,  as  discussed  earlier), 
and  are  drawn  on  the  same  set  of  axes.  Plot  symbols  will  then  be  collocated 
if  wind  and  peak  wave  directions  are  aligned. 

Because  some  events  contain  a  large  number  of  observations  of  S(f,8),  it  is 
somewhat  impractical  logistically  to  display  all  observations  for  all  events. 
Consequently,  only  certain  select  cases  are  illustrated  in  pages  following  the 
parametric  plots.  A  prestorm  condition,  all  observations  up  to  and  near  the 
peak  levels  in  a  storm,  and  several  observations  during  storm  decay  are 
displayed  for  all  events.  Shaded  symbols  in  the  graph  of  for  a  given 
event  indicate  which  observations  of  S(f,6)  are  plotted.  For  some  storms  of 
short  duration,  all  observations  are  plotted. 

Sequences  of  observed  frequency-direction  spectra,  selected  as  noted,  are 
illustrated  on  pages  following  the  parametric  plots.  For  example,  selected 
cases  from  Event  A,  indicated  by  the  shaded  symbols  on  page  A4,  are  shown 
as  sequences  of  plots  on  pages  A5  and  A6.  For  each  observation  there  are 
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four  entities  displayed:  (a)  a  three-dimensional  plot  of  the  surface  represent¬ 
ing  S(f,6),  (b)  a  linear  graph  of  the  frequency  spectrum  S(f),  (c)  a  linear 
graph  of  the  direction  spectrum  S(6)  (defined  below),  and  (d)  a  contour  plot 
of  S(f,6) .  The  same  vertical  scale  is  used  for  all  three-dimensional  plots  of 
Sif,6)  within  an  event,  so  the  relative  importance  of  each  observation  can  be 
assessed.  Prestorm  and  storm  decay  conditions  thus  show  low  relief,  and 
peak  conditions  tend  to  have  the  largest  relief. 

The  linear  plot  of  S(f)  is  shown  in  the  right  panel  of  the  three-dimensional 
box,  with  its  abscissa  parallel  to  the  frequency  axis  of  the  S(f,0)  plot.  The 
direction  spectrum  S(6)  is  the  directional  analog  of  the  frequency  spectrum 
and  shows  the  directional  distribution  of  total  energy.  Computationally,  it  is 
found  from 

S(6)  -  [A  S(f,6)df  (4) 

and  has  units  of  square  meters  per  degree.  The  analog  with  S{f)  can  be  seen 
from  the  definition  of  the  frequency  spectrum  as 

S(f)  «  f  "  S(f,6)dd  (5) 

with  units  of  square  meters  per  Hertz.  The  left  panel  of  the  three-dimensional 
box  contains  a  linear  graph  of  S(6) .  Both  S(fj  and  S(6)  are  drawn  to  constant 
scales  derived  from  the  maxima  of  these  entities  within  an  event. 

The  contour  plot  of  S(f,8)  is  shown  below  the  three-dimensional  box  for 
each  observation.  Contour  intervals  are  tenths  of  the  maximum  value  of 
S(f,0)  of  all  observations  in  an  event.  Consequently,  spectra  having  relatively 
low  energy  levels  have  both  low  relief  in  the  three-dimensional  plots  and  few, 
if  any,  lines  in  corresponding  contour  plots.  Observations  with  higher  energy 
levels  have  correspondingly  more  contour  lines. 

Because  all  four  plotted  entities  within  a  given  event  are  drawn  to  constant 
scales  within  an  event,  axis  annotation  is  rather  redundant.  It  is  also  space 
consuming.  Consequently,  axis  annotation  is  drawn  only  for  the  upper  left- 
hand  set  of  plots  on  each  page  illustrating  sequences  of  S(f,6).  Axis  annota¬ 
tion  for  all  other  plot  axes  within  an  event  is  identical. 
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4  Discussion 


The  primary  purpose  of  this  report  is  to  elucidate  the  directional  distri¬ 
butions  of  nearshore  wind  wave  energy  in  the  context  of  evolving  storms. 

With  the  illustrations  in  Appendix  A  for  guidance,  some  fundamental  proper¬ 
ties  of  storm  wave  fields  can  be  identified. 

A  primary  question  is  whether  a  storm  wave  field  can,  in  any  sense,  be 
reasonably  approximated  as  unidirectional  and  monochromatic,  or,  perhaps,  as 
corresponding  to  simple  models  of  deep-ocean  wave  fields.  A  quick  perusal 
of  Appendix  A  shows,  in  a  few  observations,  that  a  unidirectional,  monochro¬ 
matic  model  may  be  appropriate  for  some  applications  (examine,  e.g.,  Event 
R  and  conditions  near  the  peaks  of  Events  V  and  W),  but  such  conditions  are 
relatively  rare.  Furthermore,  narrow  spreads  of  naturally  occurring  spectral 
energy  in  frequency  and  direction  do  not  mean  the  waves  are  unidirectional 
and  monochromatic.  Real  ocean  waves  are  invariably  groupy  and  subject  to 
nonlinear  interactions.  These  dynamically  important  properties  can  not  be 
represented  by  unidirectional,  monochromatic  wave  models,  and  the  utility  of 
such  models  is  thereby  somewhat  restricted. 

The  wave  fields  of  Appendix  A  do  not  appear  to  conform  to  models  and 
observations  of  deep-ocean  wave  fields,  either.  For  example,  papers  by 
Mitsuyasu  et  al.  (1975)  and  Mitsuyasu  and  Uji  (1989)  indicate  that  5(f,8)  in 
the  deep  ocean  is  characterized  by  narrower  directional  spreads  near  the  spec¬ 
tral  peak  frequency  than  in  the  spectral  tails.  Spectra  shown  in  Appendix  A 
appear  generally  to  have  maximum  spreads  near  spectral  peak  frequencies. 
These  observations  mean  that  neither  a  single-wave-train  model  nor  deep- 
ocean  wave  behavior  is  appropriate  to  characterize  nearshore  storm  waves. 

While  the  spectra  of  storm  wave  fields  at  the  FRF  are  complicated,  they 
are  not  without  a  remarkable  similarity  in  structure  that  can  be  described,  at 
least  qualitatively.  In  this  part  of  the  report,  these  structural  similarities  are 
discussed,  and  some  of  the  consequent  research  questions  are  posited. 

Of  the  29  storm  events  identified  in  this  report,  25  have  S(f,6)  with  com¬ 
mon  structural  features  (described  qualitatively  below)  through  the  progress  of 
a  storm.  Four  events,  E,  P,  Q,  and  R,  do  not  follow  this  pattern.  Of  these 
four,  waves  in  Events  P  and  Q  are  associated  with  Hurricanes  Gabrielle  and 
Hugo,  respectively,  both  of  which  passed  several  hundred  kilometers  away 
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from  the  FRF.  Well -dispersed,  low-frequency  wave  trains  with  very  narrow 
directional  spreads  are  associated  with  these  events  and  are  characteristic  of 
waves  from  distant  sources.  Events  E  and  R  may  be  natural  deviates  from  the 
common  pattern  of  the  remaining  25  events. 

Storm  waves  driven  by  wind  systems  identified  in  the  FRF  annual  reports 
(Miller  et  al.  1988,  Leffler  et  al.  1989,  1990,  1991,  1992,  1993)  as  Canadian 
high-pressure  systems,  frontal  passages,  or  low-pressure  cells  forming  west  to 
south  of  the  FRF  have  frequency-direction  spectra  with  many  common  fea¬ 
tures.  Most  of  these  features  are  illustrated  in  Event  N  (pages  A47  to  A50), 
which  exemplifies  the  set  of  events.  The  progression  of  these  features  as  the 
storm  evolves  is  as  follows: 

a.  There  is  usually  a  prestorm  wave  field  consisting  of  low  energy,  long 
period  swell  with  moderate  spread  arriving  from  shore-normal  to  slight¬ 
ly  south  of  shore-normal.  For  Event  N,  this  is  the  first  case  illustrated 
on  the  top  row  of  page  A48. 

b.  During  the  primary  growth  stage,  at  the  initial  onset  of  winds,  high- 
frequency,  locally  generated  waves  appear.  Peak  period  usually  drops 
from  the  pre-storm  condition.  Peak  direction  shifts  nearly  to  align  with 
wind,  but  not  exactly  (bottom-mounted  pressure  gauges  can  detect  only 
intermediate  to  shallow  water  waves,  so  some  refraction  towards  shore- 
normal  will  take  place  even  in  short  fetches).  Peak  directional  spreads 
are  of  order  10  deg  to  20  deg,  the  smallest  of  any  time  in  the  storm. 

For  Event  N,  this  is  the  second  case  illustrated  on  the  top  row  of  page 
A48. 

c.  In  the  secondary  growth  stage  (remainder  of  the  top  row  and  first  two 
cases  of  the  second  row  of  plots  on  page  A48  for  Event  N),  peak  period 
increases,  and  the  frequency  spectrum  tends  to  conform  to  one  of  the 
conventional  nearshore  models,  such  as  TMA  (Bouws  et  al.  1985)  or 
FRF  (Miller  and  Vincent  1990).  It  is  the  directional  distribution  within 
S(f,8)  that  is  most  curious.  As  energy  at  the  peak  period  increases,  so 
does  the  characteristic  spread,  perhaps  indicating  continued  local  growth 
plus  additional  energy  radiating  in  from  growth  further  offshore.  More¬ 
over,  an  additional  lobe  of  energy  appears  at  the  peak  frequency, 
aligned  more  nearly  with  shore-normal  or,  for  Event  N,  well  south  of 
shore-normal.  Furthermore,  a  secondary  peak  of  high-frequency  energy 
(identified  as  a  ridge  in  three-dimensional  plots  of  S(f,6 )  or  a  second 
tail  off  the  main  lobe  of  energy  in  contour  plots)  appears,  making  a 
strongly  bimodal  high-frequency  tail  to  the  spectrum.  The  second 
energy  ridge  is  aligned  nearer  to  shore-normal,  and  ranges  from  20  deg 
to  40  deg  offset  from  the  initial  growth  stage  spectrum. 

d.  In  the  advanced  growth  stage  (last  two  cases  on  the  second  row  and  first 
two  cases  on  the  bottom  row  of  plots  on  page  A48  for  Event  N),  the 
spectral  peak  continues  to  move  to  lower  frequencies  and  becomes 
centered  at  very  nearly  shore-normal.  The  spread  is  rather  large,  being 
of  order  40  deg  to  50  deg,  and  the  secondary  lobe  of  low-frequency 


Chapter  4  Discussion 


energy  appears  to  be  absorbed  by  the  primary  peak.  The  structure  of 
the  bimodal,  high-frequency  tail  remains,  and  its  energy  increases. 

e.  As  storm  waves  reach  their  peak  wave  energy,  the  directionally  broad, 
energetic,  low-frequency  spectral  peak  and  directionally  bimodal,  high- 
frequency  tail  of  the  spectra  remain.  If  driving  winds  persist,  this 
structure  can  continue  relatively  unchanged  for  some  time.  In  Event  N, 
for  example,  wind  speed  remains  nearly  constant  at  16  to  18  m/sec,  and 
wind  direction  persists  at  60  deg  to  80  deg  from  shore-normal  for  more 
than  2  days  (see  page  A47).  The  sequence  of  observations  on  pages 
A48  to  A50,  showing  from  1300  on  7  Mar  89  to  1900  on  9  Mar 
89,  illustrate  a  structure  that  changes  very  little,  suggesting  that  some 
equilibrium  condition  has  been  attained. 

f.  In  the  waning  stages  of  a  storm  (the  final  plots  on  page  A50  for  spectra 
of  Event  N),  winds  usually  decrease  in  speed  and  change  direction. 
Active  local  wave  generation  abates,  and  the  high-frequency  pan  of  the 
spectrum  decays.  Low-frequency  waves  that  were  generated  offshore 
continue  to  arrive,  still  with  a  broad  directional  spread. 

g.  In  the  end  stage,  conditions  generally  return  to  the  prestorm  spectrum 
of  low-frequency,  low-energy  swell,  having  broad  directional  spread 
and  being  centered  near  the  shore-normal  direction. 

The  above  description  was  derived  specifically  from  Event  N,  but  the  other 
24  storm  events  have  the  same  features  to  a  greater  or  lesser  degree.  Even 
the  two  dissimilar  storms.  Events  E  and  R  mentioned  above,  have  some  of  the 
same  features.  The  ubiquitousness  of  these  spectral  patterns  indicates  a  high 
degree  of  consistency  in  storm  wave  evolution  at  the  FRF.  The  unexpected 
features  of  these  spectra,  low-frequency  side  lobe  of  energy  in  the  secondary 
growth  stage,  as  well  as  a  directionally  broad,  low-frequency  peak  coupled 
with  a  bimodal  high-frequency  tail  in  the  most  energetic  stages,  require  some 
explanation  so  that  guidance  appropriate  to  engineering  projects  can  be  de¬ 
rived. 

There  are  a  number  of  lines  of  approach  for  explaining  these  observations. 
Several  of  them  require  work  that  is  beyond  the  scope  of  this  report,  but  the 
nature  of  the  necessary  investigation  can  be  mentioned  here.  A  primitive 
question  is  whether  the  observed  S{f,d)  are  real  or  some  artifact  of  the  inter¬ 
action  of  the  true  wave  field  and  the  IMLE  algorithm  as  used  with  the  FRF 
linear  array  of  gauges.  If  the  wave  energy  is  all  propagating  onshore,  there  is 
no  known  reason  to  suspect  the  spectral  estimation  process.  Tests  using  syn¬ 
thetic  data  indicate  the  system  is  well  capable  of  resolving  the  broad  direction¬ 
al  spreads  and  20-deg-separated  bimodal  distributions  observed  here.  If  there 
is  wave  energy  reflected  from  nearshore,  the  directional  estimation  system 
used  in  this  report  cannot  resolve  it,  and  the  reflected  energy  will  appear  to  be 
coming  from  offshore.  This  question  can  be  addressed  with  a  more  recently 
installed  360-deg  array,  but  that  analysis  has  not  been  completed  at  this  writ¬ 
ing.  Nonetheless,  it  is  unlikely  that  reflections  will  account  for  the  present 
observations  because  the  spectral  pattern  recurs  at  many  tidal  stages  and  over 
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many  seasons,  wherein  nearshore  reflectivity  would  be  expected  to  vary  con¬ 
siderably. 

Another  possibility  is  that  linear  wave  theory  is  not  adequate  during  high- 
energy  conditions  when  the  wave  field  may  become  saturated  with  highly 
nonlinear,  actively  breaking  waves.  While  this  may  be  true  for  the  extreme 
storm  of  Event  T,  it  is  unlikely  to  be  true  for  the  very  marginally  energetic 
storm  of  Event  I,  which  was  inadvertently  included  in  this  analysis,  and  had 
an  that  exceeded  2  m  for  less  than  9  hr.  The  common  spectral  pattern 
described  above  also  occurs  in  Event  I  (see  pages  A33  and  A34),  so  it  is 
improbable  that  a  saturated  wave  field  accounts  for  these  observations. 

Two  other  possible  causes  for  these  results  are  a  common  structure  in  the 
overall  driving  wind  fields  or  the  net  effect  of  the  broad,  shallow  bathymetry 
of  the  continental  shelf  adjacent  to  the  FRF.  The  anemometer  used  in  this 
report  is  only  a  crude  indicator  of  the  larger  wind  systems  that  generate  waves 
observed  at  the  FRF.  A  possible  line  of  future  research  to  examine  the  larger 
systems  would  be  to  apply  a  wind  field  model  like  that  used 
in  the  CERC  Wave  Information  Studies  (Resio,  Vincent,  and  Corson  1982)  to 
the  storm  events  examined  in  this  report. 

In  conjunction  with  this  possible  study  would  be  the  inclusion  of  a  reason¬ 
ably  accurate  bathymetry  for  the  United  States  mid-Atlantic  continental  shelf. 
The  shelf  break,  at  depths  of  60  to  100  m,  is  80  to  100  km  distant  from  the 
FRF  (Figure  1).  Waves  with  periods  of  8  to  II  sec  (common  periods  for 
waves  observed  at  the  FRF)  begin  to  be  influenced  by  the  bottom  at  these 
depths.  Consequently,  considerable  wave  refraction  is  expected  for  intermedi¬ 
ate  and  long  waves  generated  over  the  shelf  and  in  the  deeper  Atlantic  by 
nonuniform  wind  fields  associated  with  atmospheric  fronts  and  low  pressure 
systems.  This  condition  may  account  for  the  directionally  broad,  low-frequen¬ 
cy  peaks  and  side  lobes  in  observed  storm  spectra.  However,  a  quantitative 
assessment  of  these  effects  should  be  made,  either  with  a  numerical  model  or 
with  more  extensive  observations. 

It  is  possible  that  the  directionally  bimodal,  high-frequency  spectra  may  be 
due  to  a  hummock  or  linear  bar  just  offshore  of  the  FRF  array  that  causes 
strong  refraction  for  high-frequency  waves.  This  hypothesis  can  be  tested 
with  a  detailed  bathymetric  survey  that  extends  beyond  the  area  routinely  mea¬ 
sured  by  the  FRF  staff  (generally,  the  area  shown  in  Figure  2).  Such  a  sur¬ 
vey  has  not  yet  been  done,  but  the  results  presented  herein  should  provide 
strong  motivation  for  such  work.  There  is  weak  evidence  against  a  local 
bathymetric  effect  in  the  present  data  set.  In  the  later  stages  of  Events  A,  Y, 
and  2,  winds  shifted  to  the  southeast  with  just  enough  force  to  generate  waves 
from  that  quadrant.  Though  rather  faintly  shown  in  the  figures  (see  pages  A6, 
A87,  and  A96),  bimodal  spectral  structures  appear  in  the  southeast  quadrant  in 
forms  roughly  symmetric  to  those  observed  in  the  northeast  quadrant.  It  is 
somewhat  unlikely  that  a  local  bathymetric  feature  would  induce  this  kind  of 
symmetry  in  S(f,6) . 
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Another  possible  explanation  for  the  present  observations  is  nonlinear 
interactions  among  wave  trains  during  their  generation  and  propagation  to  the 
FRF  directional  gauge.  Various  interaction  mechanisms  have  been  proposed 
by  the  WAMDI  Group  (1988),  Freilich,  Guza,  and  Elgar  (1990),  Resio  and 
Perrie  (1991),  and  others.  Testing  these  hypotheses  requires  more  detailed 
spatial  measurements  of  S(f,6 )  than  are  currently  available.  However,  an 
experiment  can  and  should  be  designed  for  this  purpose  if  the  current  results 
are  to  be  fully  explained. 

By  way  of  guidance  for  model  testing  and  engineering  design,  a  representa¬ 
tive  example  of  S(f,6)  for  late  stage  growth  and  peak  storm  conditions  is 
shown  in  Figure  4.  It  is  simply  an  average  of  the  19  observations  from  the 
relatively  stationary  peak  storm  conditions  of  Event  N.  It  contains  the  prima¬ 
ry  features  noted  above  and  represents  a  nearly  stationary  dynamic  process 
that  continued  for  56  hr.  It  should  be  noted  that  although  the  high-frequency 
tail  has  a  spectral  density  less  than  the  spectral  peak,  the  total  energy  in  the 
part  of  the  tail  that  is  bimodal  (found  from  the  area  under  the  frequency  spec¬ 
trum  over  the  high-frequency  range  of  frequencies)  is  roughly  25  percent  to 
50  percent  of  the  total  energy.  Because  the  bimodal  tail  occurs  at  peak  energy 
conditions,  there  is  considerable  energy  in  the  tail  region. 

The  key  point  is  that  a  wave  generation  and  transformation  model  at  conti¬ 
nental  shelf  and  larger  scales  clearly  must  have  a  result  that  looks  qualitatively 
like  Figure  4  to  represent  conditions  at  the  FRF  with  reasonable  verity.  Addi¬ 
tionally,  nearshore  process  models  of  the  FRF  will  be  improved  if  a  spectrum 
like  that  in  Figure  4  is  used  to  describe  the  wave  boundary  condition  outside 
the  surf  zone  in  high-energy  cases. 
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Mean  Frequency-Direction  Spectrum 
19  Cases  from  1300,  7  Mar  89  to  1900,  9  Mar  89 
Hmo  “  3.72  ±  0.20  m 
depth  *  8.74  ±  0.50  m 
Tp  —  11.98  sec  8 p  *=  4.0  deg 


Figure  4.  A  characteristic  storm  spectrum  from  the  FRF 
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This  report  elucidates  the  nature  of  frequency-direction  spectra  evolving 
during  storm  events  occurring  at  the  FRF.  The  primary  result  is  a  series  of 
graphic  presentations  of  spectra  measured  with  a  high-resolution  directional 
wave  gauge  in  an  8-m  depth  approximately  900  m  offshore.  Data  from  a 
5-year  observation  period  have  been  searched  for  storm  events  defined  Kv  a 
characteristic  wave  height  in  excess  of  2  m  at  two  points  at  least  24  hr  -art 
in  the  same  storm  system.  A  total  of  29  events  have  been  isolated  for  closer 
scrutiny.  Time  series  plots  of  climatic  parameters  and  sequences  of  plots  of 
S(f,8)  are  used  to  reveal  the  evolutionary  nature  of  storm  wave  fields. 

The  data  indicate  that  there  are  two  basic  types  of  events:  (a)  high  waves 
from  strong  but  distant  sources,  and  (b)  localized  wave  generation  from  Cana¬ 
dian  high-pressure  systems,  frontal  passages,  and  atmospheric  low-pressure 
cells  passing  the  FRF  from  the  west  or  south.  Of  the  29  events  identified,  27 
are  of  this  latter  type.  Of  the  27  localized  events,  25  demonstrated  very  com¬ 
mon  patterns  of  spectral  evolution.  The  patterns  include  initial  localized 
waves  aligned  nearly  with  the  local  wind,  directional  broadening  of  spectral 
peaks  as  the  wave  fields  grow,  and  a  bifurcation  of  high-frequency  spectral 
tails  into  two  modes:  (a)  one  nearly  aligned  with  the  wind,  and  (b)  another, 

20  deg  to  40  deg  away,  aligned  nearer  to  the  shore-normal  direction.  One- 
quarter  to  one-half  of  the  total  energy  is  contained  in  the  bimodal  tails,  a 
consideration  of  consequence  because  such  spectra  occur  at  the  peaks  of 
storms,  and  because  the  beating  of  such  spectral  waves  against  the  shoreface 
can  affect  such  nearshore  processes  as  the  patterns  of  beach  erosion  and  infra- 
gravity  wave  generation. 

The  main  conclusions  here  are  that  storm  waves  occur  at  the  FRF  with 
curious  spectral  shapes,  and  that  accounting  for  the  processes  that  give  rise  to 
those  shapes  is  important.  It  is  suggested  that  a  more  detailed  coverage  of 
wind,  wave,  and  bathymetric  observations  over  100-km  scales  from  the  FRF 
is  warranted.  It  is  certain  that  nearshore  processes  driven  by  waves  having 
such  elaborately  structured  spectra  will  be  different  from  estimates  generated 
by  over-simplified  wave  models  (typically  unidirectional  and  monochromatic) 
as  even  now  still  haunts  engineering  guidance.  There  is  good  reason  to  be¬ 
lieve  the  spectra  reported  here  are  real  and  not  some  artifact  of  the  analysis 
method  used.  Use  of  these  results  for  testing  wave  generation  and  transforma¬ 
tion  models  at  scales  of  the  continental  shelf  width,  and  as  seaward  wave 
boundary  conditions  for  nearshore  process  models,  will  certainly  yield  more 
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realistic,  meaningful,  and  accurate  representations  of  nature,  and  thereby  an 
improved  coastal  engineering  science. 
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wind  speed  {m/sec)  direction  (deg)  A0  (deg)  depth  (m)  T  (sec)  H  (m) 

0  18  -180  0  180  0  40  7.0  l.l  0  12  0.0  3.2 


Climatological  Parameters 
Event  A:  9  Oct  86  to  14  Oct  86 
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wind  speed  (m/sec)  direction  (deg)  A 0  (deg)  depth  (m)  T  (sec) 

0  13  -ISO  0  180  0  46  7.1  8.8  0  * 


Climatological  Parameters 
Event  B:  17  Oct  86  to  21  Oct  86 
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wind  speed  (m/sec)  direction  (deg)  A 6  (deg)  depth  (m)  T  (sec) 

0  18  -ISO  0  ISO  O  80  7.2  S.4  0  * 


Climatological  Parameters 
Event  C:  15  Feb  87  to  20  Feb  87 


Appendix  A  Event  Time  Series  end  Spectral  Sequences 


Appendix  A  Event  Time  Series  end  Spectrel  Sequences 


All 


depth  (m) 


Climatological  Parameters 
Event  D:  9  Mar  87  to  18  Mar  87 
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wind  speed  (m/sec)  direction  (deg)  t0  (deg)  depth  (m)  T  (sec) 

o  11  -ISO  0  180  0  55  7.4  8.8  0  P 


Climatological  Parameters 
Event  E:  22  Mar  87  to  26  Mar  87 
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0  18  -180  0  180  0  46  7.t  8.6  0 


Climatological  Parameters 
Event  3:  11  Oct  87  to  17  Oct  87 
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wind  speed  (m/$ec)  direction  (deg)  A 8  (deg)  depth  (m)  T  (sec! 

O  12  >180  0  180  0  60  7.0  8.4  0 


Climatological  Parameters 
Event  G:  26  Nov  87  to  1  Dec  87 
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wind  speed  (m/sec)  direction  (deg)  A 6  (deg)  depth  (m)  T  (sec) 

0  13  -ISO  0  ISO  0  45  7.3  8.4  0  P 


Climatological  Parameters 
Event  H:  14  Dec  88  to  17  Dec  88 
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wind  speed  (m/sec)  direction  (deg)  A 8  (deg)  depth  (m)  T  (sec) 

0  13  -ISO  0  ISO  0  45  7.3  B.S  0 


Climatological  Parameters 
Event  I:  3  Jan  89  to  6  Jan  89 
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wind  speed  (m/sec)  direction  (deg)  A 8  (deg)  depth  (m)  T  (sec) 

0  19  >190  0  ISO  0  66  7.0  MO  ’ 


Climatological  Parameters 
Event  L:  22  Feb  89  to  26  Feb  89 
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wind  speed  (m/sec)  direction  (deg)  A 6  (deg)  depth  (m)  T  (sec) 

0  15  -ISO  0  ISO  0  65  7  2  8  6  0  " 


Climatological  Parameters 
Event  0:  21  Mar  89  to  25  Mar  89 
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wind  speed  (m/sec)  direction  (deg)  hB  (deg)  depth  (m)  T  (sec 

0  17  -180  0  ISO  0  36  7.2  8  8  0  * 


Climatological  Parameters 
Event  Q:  18  Sep  89  to  28  Sep  89 
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wind  speed  (m/sec)  direction  (deg)  LB  (deg)  depth  (m)  T  (sec) 

0  18  -180  0  180  0  60  7.2  8.8  0  P 


Climatological  Parameters 
Event  U:  6  Mar  90  to  9  Mar  90 


Appendix  A  Event  Time  Series  end  Spectral  Sequences 


A73 


Appendix  A  Event  Time  Series  end  Spectrel  Sequences 


A75 


Appendix  A  Event 


Frequency-Direction  Spectral  Plots:  Event  W 


0400.  17  Nov  90 
Hmo  —  0.22  m 


1300,  17  Nov  90 
Hmo  —  131  m 


1600.  17  Nov  90 
Hmo  —  146  m 


1900.  17  Nov  90 
Hmo  ”  1.79  m 


wmmmm 


2200.  17  Nov  90 
Hmo  —  2.45  m 


0100.  18  Nov  90 
Hmo  ■*  2.51  m 


0400.  18  Nov  90 
Hmo  “  2.39  m 


0700.  18  Nov  90 
Hmo  “  2.44  m 


ip**. 


afim 


1000.  18  Nov  90 
Hmo  <■  2.42  m 


1300,  18  Nov  90 
Hmo  —  2.13  m 


1900.  18  Nov  90 
Hmo  “  1.99  m 


0100,  19  Nov  90 
Hmo  -  2.04  m 


.jRC- 


m 


-7  ■  * 


Appendix  A  Event  T*ne  Series  and  Spectral  Sequences 


wind  speed  (m/sec)  direction  (deg)  A0  (deg)  depth  <m)  T  (sec)  H  (m) 

0  12  -ISO  0  180  0  40  7.1  8.6  0  14  0.0  2.8 


Climatological  Parameters 
Event  X:  7  Dec  90  to  10  Dec  90 
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wind  speed  (m/sec)  direction  (deg)  A 8  (deg)  depth  (m)  T  (sec) 

0  16  -180  0  ISO  0  38  7.3  1.6  0  " 


Climatological  Parameters 
Event  Z:  18  Apr  91  to  21  Apr  91 
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wind  speed  (m/sec)  direction  (deg)  A 8  (deg)  depth  (m)  T  (sec) 

0  14  -ISO  0  J BO  O  3 B  7.4  boo  p 


Climatological  Parameters 
Event  2:  23  Jun  91  to  28  Jun  91 
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d  Water  depth 

/  Cyclic  frequency 

fK  Cyclic  frequency  of  n*  component  wave 

f  Spectral  peak  frequency 

fN  High-frequency  boundary  for  a  discrete  spectrum 

fx  Low-frequency  boundary  for  a  discrete  spectrum 

H ^  Spectrum-based  characteristic  wave  height 

I(f,d )  Direction  integral  of  directional  distribution  function 

IMLE  Iterative  Maximum  Likelihood  Estimation 

m  Index  integer 

M  Upper  limit  of  indices  denoted  m 

MLE  Maximum  Likelihood  Estimation 

n  Index  integer 

N  Upper  limit  of  indices  denoted  n 

S  Spectral  density 

S(f)  Frequency  spectral  density 

S(fJ  Frequency  spectral  density  at  n*  discrete  frequency 
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S(0) 
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Frequency-direction  spectra)  density 

Frequency-direction  spectral  density  at  n*  discrete 
frequency  and  m*  discrete  direction 

Direction  spectral  density 

Characteristic  wave  period  derived  from  spectral  peak 

Convergence  exponent  in  IMLE  method 

Convergence  coefficient  in  IMLE  method 

Directional  spread  of  wave  energy  at  spectral  peak 
frequency 

Finite  frequency  increment 
Wave  direction 

Element  m  of  a  discrete  set  of  wave  directions 
Upper  bound  of  discrete  direction  coordinate 
Peak  wave  direction 

Lower  bound  of  discrete  direction  coordinate 
Angle  subtending  25  percent  of  wave  energy 
Angle  subtending  75  percent  of  wave  energy 
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